Formation of Cooper pairs in quantum oscillations of electrons in plasma 
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We study low energy quantum oscillations of electron gas in plasma. It is shown that two electrons 
. . , participating in these oscillations acquire additional negative energy when they interact by means 

of a virtual plasmon. The additional energy leads to the formation a Cooper pair and possible 
»' I existence of the superconducting phase in the system. We suggest that this mechanism supports 

^— »' , slowly damping oscillations of electrons without any energy supply. Basing on our model we put 

forward the hypothesis the superconductivity can occur in a low energy ball lightning. 
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INTRODUCTION 



In the majority of cases plasma contains a great deal of electrically charged particles. Free electrons in metals can be 
treated as a degenerate plasma. If the motion of charged particles in plasma is organized, one can expect the existence 
^ ■ of macroscopic quantum effects. For instance, collective oscillations of the crystal lattice ions can be exited in metals. 
These oscillations are interpreted as quasi-particles or phonons. Phonons are known to be bosons. The exchange of 
a phonon causes the effective attraction between two free electrons in plasma of metals. This phenomenon is known 



^ as the formation of a Cooper pair ^] and is important in the explanation of the superconductivity of metals. Cooper 



pairs can exist during a macroscopic time resulting in an undamped electric current. Unfortunately the formation 



^ ' of Cooper pairs is possible only at very low temperatures. Nowadays there are a great deal of attempts to obtain a 



superconducting material at high temperatures 



^ • Without an external field, e.g. a magnetic field, the motion of electrons in plasma is stochastic. It results in rather 
short times of plasma recombination. The laboratory plasma at atmospheric pressure without an energy supply 
recombines during ^ 10~^s @. That is why free oscillations of electrons in plasma, i.e. without any external source 
of energy, will attenuate rather fast. Therefore to form a Cooper pair in plasma one should create a highly structured 
electrons motion to overcome the thermal fluctuations of background charged particles. Nevertheless it is known that 
a plasmon superconductivity can appear [J. 

In Ref. [5*] we studied quantum oscillations of electron gas in plasma on the basis of the solutions to the non-linear 
Schrodinger equation. The spherically and axially symmetrical solutions to the Schrodinger equation were found in 
our works. We revealed that the densities of both electrons and positively charged ions have a tendency to increase in 
the geometrical center of the system. The found solutions belong to two types, low and high energy ones, depending 
^\ ^ on the branch in the dispersion relation (see also Sec. [TT|- We put forward a hypothesis that a microdose nuclear 
■ fusion reactions can happen inside a high energy oscillating electron gas. This kind of reaction can provide the energy 
^ I supply to feed the system. In our works we suggested that the obtained solutions can serve as a model of a high 
• I— I . energy ball lightning. Note that the suggestion that nuclear reactions can take place inside a ball lightning was also 
^ ' put forward in Ref. . 

^ , In the present work we study a low energy solution (see Sec. |TT| which is also presented in our model As we 
- - ' mentioned above a non-structured plasma usually recombines during several milliseconds. It is the main difficulty 
in constructing of long-lived plasma structures. In a high energy solution internal nuclear fusion reactions could 
compensate the energy losses. We suggest that the mechanism, which prevents the decay of low energy oscillations 
of the electron gas, is based on the possible existence of the plasma superconductivity. In frames of our model we 
show that the interaction of two electrons, described by spherical waves, can be mediated by a virtual plasmon. For 
plasmon frequencies corresponding to the microwave range, this interaction results in the additional negative energy. 
Then it is found that for a low energy solution with the specific characteristic this negative energy results in the 
effective attraction. Therefore a Cooper pair of two electrons can be formed because the interacting electrons should 
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have oppositely directed spins. Using this result it is possible to state that the friction occurring when electrons 
propagate through the background matter can be significantly reduced due to the superconductivity phenomenon. 
We also discuss the possible applications of our results for the description of the a low energy ball lightning. It should 
be noticed that the appearance of the superconductivity inside of a ball lightning was considered in Ref. [7[ . 

This paper is organized as follows. In Sec. [Til we briefly outline our model of quantum oscillations of the electron 
gas which was elaborated in details in Ref. [5]. Then in Sec. IIIII we study the interaction of two electrons when 
they exchange a virtual plasmon and examine the conditions of the effective attraction appearance. In Sec. IIVI we 
summarize and discuss the obtained results. 



II. BRIEF DESCRIPTION OF THE MODEL 



In Ref. Q we discussed the motion of an electron with the mass m and the electric charge e interacting with 
both background electrons and positively charged ions. The evolution of this system is described by the non-linear 
Schrodinger equation, 



where 



dilj 



(2.1) 



(2.2) 



is the term which is responsible for the interaction of an electron with background matter which includes electrons 
and positively charged ions with the number density ni{r,t). In Eq. (|2.ip the wave function is normalized on the 
number density of electrons, [■(/'(r, t)p = 7ie(r,t). 

Note that in Eq. (|2.ip the c-number wave function ip = ipi^jt) depends on the coordinates of a single electron. 
Generally in a many-body system this kind of approximation is valid for an operator wave function, i.e. when it is 
expressed in terms of creation or annihilation operators. In Ref. it was shown that it is possible to rewrite an 
exact operator Schrodinger equation using c-number wave functions for a single particle. Of course, in this case a set 
of additional terms appears in a Hamiltonian. Usually we neglect these terms (see Ref. Q). 

The solution to Eqs. (|2.ip and (|2.2p can be presented in the following way: 



(2.3) 



where hq is the density on the rim of system, i.e. at rather remote distance from the center. The perturbative wave 
fimction (5V'(r,f) for the spherically symmetrical oscillations has the form (see Ref. Q), 



sin kr 



(2.4) 



where Ak is the normalization coefficient. 

The dispersion relation which couples the frequency of the oscillations uj and the quantum number k reads 



4 = 



1± 1-4-^^ 



(2.5) 



It is convenient to rewrite Eq. (|2.5|) in the new variables x = Lu/u>p and y — k^Jh/2ujpm. Now it is transformed into 
the form. 



y± 



1± 1 



1/2- 



(2.6) 



We represent the functions y± versus x on Fig. [TJ One can see on Fig. [T] that there are two branches y± in the 
dispersion relation, upper and lower ones. It is possible to identify the upper branch with a high energy solution and 
the lower branch - with a low energy one. 

Generally Eqs. (|2.ip and (|2.2p include the interaction of the considered electron with all other electrons in the 
system. However since we look for the solution in the pertubative form (|2.3p and linearize these equations, the 
solution to Eqs. (12. 4p and (|2.5p describe the interactions of the considered electron only with background electrons, 
which correspond to the coordinate independent wave function ip^, rather than with electrons which participate in 
spherically symmetrical oscillations. 
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FIG. 1: The dispersion relation for quantum oscillations of electron gas in plasma based on the results of Ref. [^. 

III. INTERACTION OF TWO ELECTRONS PERFORMING OSCILLATIONS CORRESPONDING TO 

THE LOW ENERGY BRANCH 



In this section we discuss the interaction between two electrons participating in spherically symmetrical oscillations. 
This interaction is supposed to be mediated by a plasmon field <p. To account for the electron interaction with a 
plasmon one has to include the additional term eip into the Hamiltonian in Eq. (|2.ip . We also admit that electron 
gas oscillations correspond to the lower branch y_ in the dispersion relation (|2.6p . i.e. we are interested in the 
dynamics of a low energy solution. 

A plasmon field distribution can be obtained directly from the Maxwell equation for the electric displacement field 
D, 



(VD) = 47rp, 



(3.1) 



where p is the free charge density of electrons. As it was mentioned in Ref. [H*] the vector potential for spherically 



symmetrical oscillations can be put to zero A 
potential 95 in the form of the Fourier integral, 

we get the formal solution to Eq. (j3.ip as 



0. Therefore only longitudinal plasmons can exist. Presenting the 



dcjd^p 



</j(p,a;)e" 



47r y df'd3r'G(r - 



where 



G(r,t) 



dwd^p 



(2^)4 p2£,(p,c^) 



,t-i')p(r',i'), 



^— icjt+ipr 



(3.2) 



(3.3) 



(3.4) 



is the Green function for the longitudinal plasmon. In Eq. ()3.4p £1 stands for the longitudinal permittivity. Note 
that to choose the exact form of the Green function (|3.4p one has to specify passing-by a residue at p = 0. We will 
regularize the Green function by putting a cut-off at small momenta. It should be noticed that the more rigorous 
derivation of the plasmon Green function is presented in Ref. Q . 

To relate Eq. (|3.3p with the distribution of electrons we can take that free charge density is proportional to the 
probability distribution, p = e\4'\'^. Then, using Eq. (|3.3p and the Schrodinger equation (j2.1D with the new Hamiltonian 
H = Hq + eif we obtain the closed non-linear equation for the the electron wave function only. 
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Now we write down the new Hamiltonian explicitly, 

H = Ho + V, V ^iire^ J dt'd^r'G{r~r',t^t')\'iP{r',t')\^. (3.5) 

As we mentioned in Sec. |TT] the term Hq describes the properties of a single electron including its interactions with 
background particles. The second term V in the Hamiltonian p.Sp is responsible for the interactions between two 
separated electrons. The scattering of two electrons is due to this interaction. 

It is possible to introduce the additional energy acquired in scattering of two electrons as 

AE= ^ine^ J dt'd^r'dtd^r^* {r,t)i:* {r' ,t')G{r ~ r' ,t ~ t')i:{r,t)il;{r' ,t'). (3.6) 

In Eq. (|3.6p we add one additional integration over t which is convenient for further calculations. To take into account 
the exchange effects one has to replace the products of the electrons wave functions in the integrand in Eq. (|3.6p with 
the symmetric or antisymmetric combinations of the basis wave functions, 

V'(r, t') ^-L[e--i*--^*>,^ (r)^'^, (r') ± e--^*--^*^,., {v)4>u, (r')], 

r{v,t)r{v\t') -.-L[ei/i*+i^*>;(r)V;(r') ± e'^^^+'^^'V; (r)^;^ (r')], (3.7) 

where fci^2 and qi.2 are the initial and final quantum numbers, and (^1.2 and /1.2 are the initial and final frequencies 
of electron gas oscillations. Since the total wave function of two electrons, which also includes their spins, should be 
anti-symmetric, the coordinate wave functions in Eq. (|3.7p which are symmetrical correspond to anti-parallel spins. 
For the same reasons the anti-symmetrical functions in Eq. p.7p describe particles with parallel spin. 

It is worth to mention that in calculating the additional scattering energy p.6p electrons should be associated 
with the perturbative wave function (|2.4p rather than with the total wave function (|2.3p . However we will keep 
the notation ip for an electron wave function to avoid cumbersome formulas. Using the basis wave functions (|2.4p . 
Eqs. p.6p and p.7p as well as the known value of the integral. 



[°° ds IT 

I — sin(ps) svaiks) sin(<7s) = — [sign(p k — q) + sign(p + g — fc) — sign(p k + q) — sign(p ~ k ~ q)], (3.8) 
Ja s 8 

we obtain for the additional energy the following expression: 

e^TT^ /-^ dp 



AE -<5(/i + /2 - - iu,)A;^A;^Ak,Ak, ^ ^ 





[sign(p + fci - gi) + sign(p + gfi - fci) 



- sign(p -I- fci -I- 91) - sign(p - fci - qi)] x [1 ^ 2] 

± —, — J- r [sign(p + A:2 - gi) + sign(p + qi ~ fc2) 

- sign(p + fc2 + qi) - sign(p - k2 - qi)] x [1 ^ 2] I, (3.9) 



where the symbol [1 <-> 2] stands for the terms similar to the terms preceding each of them but with all quantities 
with a subscript 1 replaced with the corresponding quantities with a subscript 2. 

We suggest that electrons move towards each other before the collision and in the opposite directions after the 
collision, i.e. ki = —k2 = k and qi — —(72 = q- In this situation the electrons occupy the lowest energy state (see, 
e.g., Ref. [l3|)- Taking into account that the frequencies and the coefficients Ak are even functions of k [see Eq. ()2.5p ] 
we rewrite the expression for the additional energy in the form, 



2 4 roc 

AE =Sif-u;)\A,f\Ak\'^ 



dp 



16 Jq p'^ ei{p, f - uj) 
X [sign(p + k- q) + sign(p + q - k) - sign(p + k + q) - sign(p - fc - q)f, (3.10) 

where uj — uji ~ lu2 and / = /i = /2- It is important to note that only symmetric wave functions from Eq. p.7p 
contribute to Eq. p.lOp . The anti-symmetric wave functions convert Eq. (|3.9p to zero. It means that only electrons 
with oppositely directed spins contribute to the additional energy. 
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Now we should eliminate the energy conservation delta function. One can make it in a standard way, 5{f — uj) — > 
Time/(27r), and then divide the whole expression on the observation time. Since the energy is concerned in the 
scattering one has two possibilities for the momenta, k + q ~ {) oi k — q = Q. We can choose, e.g., the latter case 
supposing that A; > 0. In this situation one has for Ai? the following expression: 

2 3 /-oo 1 -, 

^^[l-sign(p-2fc)P. (3.11) 
Finally accounting for the step function in the integrand we can present the additional energy in the form, 

^E^\A.f2^r%^. (3.12, 



P^£i(P,0)' 

Depending on the sign of the the additional energy one has either effective repulsion or effective attraction. 

As we have found in Ref. [5i] the typical frequencies of electron gas oscillations in our problem ~ 2ajp, which is 
~ 10"'^'^ s~^. Such frequencies lie deep inside the microwave region. The permittivity of plasma for this kind of 
situation has the following form (see, e.g., Ref. [Tlj'): 

eK^) = l-^^, (3.13) 

where v is the transport collision frequency. In Eq. (|3.13p we neglect the spatial dispersion. It is clear that in Eq. p.l2p 
one has to study the static limit, i.e. w ^ 0. Supposing that the density of electrons is relatively high, i.e. uop 3> 
one gets that the permittivity of plasma becomes negative and the interaction of two electrons results in the effective 
attraction. 

Note that at the absence of spatial dispersion the integral is divergent on the lower limit. Such a divergence 
is analogous to the infrared divergence in quantum electrodynamics. We have to regularize it by the substitution 
i/cLD, where an — yAirehi^Jki^ is the Debye length and T is the plasma temperature. 

Besides the energy A£^, which is acquired in a collision, an electron also has its kinetic energy and the energy of 
the interaction with background charged particles. For the effective attraction to take place we should compare the 
additional negative energy with the eigenvalue of the operator iJp, which includes both kinetic term and interaction 
with other electrons and ions. We suggest that two electrons, which are involved in the scattering, move towards each 
other, i.e. lui = u)2 = io. Therefore, if we get that 2fujj + Ai? < 0, the effective attraction takes place. 

Calculating the integral in Eq. (|3.12p one obtains the following inequality to be satisfied for the attraction to happen: 



^^TT^ nov^hv [y — y) ^ 1 I hu!. 

2/5 ' ^"2V2^u 



where v = \J ksT jm. is the thermal velocity of background electrons in plasma and ks ~ 1.38 x 10~^^erg/K is the 
Boltzmann constant. In Eq. p.l4p we use the variables x and y which were introduced in Sec. [Tll We also assume 
that = ric/k'^, where ric is the density in the central region [see Eq. (|2.4p ]. As it was predicted in our work 0], 
ric should have bigger values that the density on the rim of the system. We can also introduce two dimensionless 
parameters, ^ and a, to relate the density in the center, the density of free electrons, and the rim density: ric = ^?io 
and Uf. — ariQ. 

We suppose that a ^ 0.1. It is known (see, e.g., Ref. [i3|) that the number density of free electrons in a spark 
discharge is equal to 10^^ — 10^^ cm^'^. Approximately the same electrons number densities are acquired in a streak 
lightning. If one takes the rim density equal to the Loschmidt constant no — 2.7 x 10^^ cm~^, we arrive to the maximal 
value of a ~ 0.1. In this case mainly the collisions with neutral atoms would contribute to the transport collision 
frequency. One can approximate this quantity ds v ^ varic, where a = ttR^ is the characteristic collision cross section 
and R ~ 10~® cm is the typical atomic radius. Note that the process in question happens in the center of the system. 
That is why we take the central number density Uc in the expression for v. 

Finally we receive the following inequality to be satisfied for the effective attraction between two electrons to take 
place: 



^ w 1.1 X 10 , 3.15) 

256 mujQ 
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FIG. 2: The allowed regions in the x) plane for the existence of the effective attraction between two electrons performing 
quantum oscillations corresponding to the lower branch on Fig. [1] We also add one more vertical axis on the right hand side 
demonstrating the behaviour of the quantum number k. Three zones for the temperatures T = 500 K, 600 K and 700 K are 
shown here. 

where r = T/Tq. In Eq. (|3.15p we take the reference temperature To = 700 K. Now one calculates other reference 
quantities like fo and loq as vq = ^Jkj^I\^lm ~ 10^ cm/s and loq — Ane^nQ/m k, 2.9 x lO^"* s~^. 

As one can see on Fig. [1] [see also Eq. (|2.6p ] the parameter y is always more than 1 for the upper branch. It means 
that the inequality in Eq. (j3.15p can never be satisfied. On the contrary, there is a possibility to implement the 
condition (|3.15|1 for the lower branch on Fig. [TJ It signifies the existence of the effective attraction between electrons. 

On Fig. [2] we plot the allowed regions where the effective attraction happens. This plot was created for the constant 
value of a = 0.1. We analyze the condition (|3.15p numerically. The calculation is terminated when we reach the point 
x = 2. It has the physical sense since oscillations cannot exist for x <2 according to Eq. (|2.6p (see also Fig. [l}. 

One can notice that always there is a critical value of the central density. It is equal to 120no, 105no and 95no 
for the temperatures T = 500 K, 600 K and 700 K respectively. It means that in the center of the system there is an 
additional pressure of about 100 atm. 

One has to verify the validity of the used approach checking that y > y in Eq. (|3.14p . It would signify that that 
the the quantum number k is more than l/2ao. The values of y are 0.42, 0.38 and 0.35 for T = 500 K, 600 K and 
700 K respectively. Comparing these values with the curves on Fig. [5] one can see that the condition y > y is always 
satisfied for all the allowed zone of the effective attraction existence. 

For the effective attraction to happen one should have the negative values of £; (0) in Eq. (|3.12p . Using Eq. p.l3p 
one can conclude that the following inequality should be satisfied: v < ujp or equivalently ^ < 3.3 x lO^^/aJr. For 
a — 0.1 and temperatures T = 500 K, 600 K and 700 K we get the critical values of £, as 1.2 x 10^, 1.1 x 10^ and 
1.0 X 10'^ respectively. One sees on Fig. [5] that all the curves lie far below the critical values of Therefore the 
condition < is always satisfied. 

IV. SUMMARY AND DISCUSSIONS 

In conclusion we mention that we studied spherically symmetrical quantum oscillations of electron gas in plasma 
which correspond to a low energy solution to the Schrodinger equation [g]. We suggested that two independent 
plasma excitations, corresponding to two separate electrons, can interact via exchange of a plasmon. Since the typical 
oscillations frequencies are very high and lie in the deep microwave region the plasma permittivity is negative and 
the interaction of two electrons result in the appearance of an attraction. For the this attraction to be actual we 
compared it with the kinetic energy of electrons and the energy of interaction of electrons with other background 



7 



charged particles. The total energy of an electrons pair turned out to be negative for the lower branch in the dispersion 
relation (see Fig. [T] and Ref. [5|]). Although we did not demonstrate it explicitly, using the same technique as in Sec. lIIII 
one can check that the total energy of electrons is always positive for the upper branch. It means that even additional 
negative energy does not cause the attraction of electrons. It is interesting to mention that to obtain the master 
equations (|2.ip and (|2.2p in Ref. [5| we neglected the exchange effects between an electron which corresponds to the 
wave function ip and the rest of background electrons (see also Ref. [1]). This crude approximation is valid if the 
superconductivity takes place because in this case the friction, i.e. the interaction, between oscillating and background 
electrons is negligible. 

The proposed plasmon superconductivity happening at spherically symmetrical oscillations of electrons could be 
implemented in a low energy ball lightning. Theoretical and experimental studies of ball lightnings have many years 
history (see, e.g., Ref. [TsjV This natural phenomenon, happening mainly during a thunderstorm, is very rare and we 
do not have so many its reliable witnesses. There were numerous attempts to generate in a laboratory stable structures 
similar to a ball lightning. Many theoretical models aiming to describe the observational data were put forward. In 
putting forward a ball lightning model one should try to explain these properties on the basis of the existing physical 
ideas without involvement of extraordinary concepts (see, e.g. Ref. [15']) though they look very exciting. However 
none of the available models could explain all the specific properties of a fireball. 

Among the existing ball lightning models one can mention the aerogel model [3]. According to this model fractal 
fibers of the aerogel can form a knot representing the skeleton of a ball lightning. Using this model it is possible 
to explain some of the ball lightning features. The interesting model of a ball lightning having complex onion-like 
structure with multiple different layers was put forward in Ref. [l^. However in this model an external source of 
electrical energy is necessary for a ball lightning to exist for a long time. The hypothesis that a ball lightning is 
composed of molecular clusters was described in Ref. [l^ • This model can explain the existence of a low energy ball 
lightning. In Ref. [18] the sophisticated ball lightning model was proposed, which is based on the non-trivial magnetic 
field configuration in the form of closed magnetic lines forming a knot. In frames of this model one can account for the 
relatively long life time of a ball lightning. It is worth to be noticed that though the properties of low energy fireballs 
can be accounted for within mentioned above models they are unlikely to explain its regular geometrical shape. The 
author of the present work has never observed a ball lightning, however according to the witnesses a fireball resembles 
a regular sphere. There are many other ball lightning models which are outlined in Ref. (Toj . 

The most interesting ball lightning properties are (see, e.g., Ref. 

• There are both high and low energy ball lightnings. The energy of a low energy ball lightning can be below 
several hundred kJ. However the energy of a high energy one can be up to 1 MJ. 

• There are witnesses that a ball lightning can penetrate through a window glass. Sometimes it uses existing 
microscopic holes without destruction of the form of a ball lightning. It means that the actual size of a ball 
lightning is rather small and its visible size of several centimeters is caused by some secondary effects. Quite 
often a ball lightning can burn tiny holes inside materials like glass to pass through them. It signifies that the 
internal temperature and pressure are very high in the central region. 

• A fireball is able to follow the electric field lines. It is the indirect indication that a ball lightning has electro- 
magnetic nature, e.g., consists of plasma. 

• A ball lightning has very long life-time, about several minutes. In Sec.|I]we mentioned that unstructured plasma 
without external energy source looses its energy and recombines extremely fast. If we rely on the plasma models 
of a ball lightning, it means that either there is an energy source inside of the system or plasma is structured 
and there exists a mechanism which prevents the friction in the electrons motion. 

• In some cases there were reports that a ball lightning could produce rather strong electromagnetic radiation 
even in the X-ray range. It might signify that energetic processes, e.g. nuclear fusion reactions, can happen 
inside a fireball. 

Quantum oscillations of electron gas in plasma 5] can be a suitable model for a ball lightning. For example, it 
explains the existence of two types of ball lightnings, low and high energy ones. Our model accounts for the very 
small ~ 10~^ cm and dense central region as well as indirectly points out on the possible microdose nuclear fusion 
reactions inside of a high energy ball lightning. 

In the present paper we suggested that superconductivity could support long life-time of a low energy fireball and 
described that this phenomenon could exist in frames of our model. It is unlikely that nuclear fusion reactions can 
take place inside a low energy ball lightning, i.e. this type of a ball lightning cannot have an internal energy source. 
Moreover as we mentioned in Ref. [5|, electrons participating in spherically symmetrical oscillations cannot emit 
electromagnetic radiation. It means that the system does not loose the energy in the form of radiation. As a direct 
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consequence of our model we get the regular spherical form of a fireball gratis. In Sec. IIIII w e made the numerical 
simulations for the temperature range 500 — 700 K. According to the fireball witnesses [ij] some low energy ball 
lightnings do not combust materials like paper, wood etc. The combustion temperature of these materials lie in the 
mentioned above range. It justifies our assumption. 

Note the a ball lightning seems to be a many-sided phenomenon and we do not claim that our model explains all the 
existing electrical atmospheric phenomena which look like a ball lightning. However in our opinion a certain class of 
fireballs with the above listed properties is satisfactory accounted for within the model based on quantum oscillations 
of electron gas in plasma. 
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